Chaos or Order in Double Barred Galaxies? 



Witold Maciejewski 

INAF - Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy, 
and Obserwatorium Astronomiczne Uniwersytetu Jagiellonskiego, Krakow, Poland 



Abstract. Bars in galaxies are mainly supported by particles trapped around closed 
periodic orbits. These orbits respond to the bar's forcing frequency only and lack free 
oscillations. We show that a similar situation takes place in double bars: particles get 
trapped around orbits which only respond to the forcing from the two bars and lack 
free oscillations. We find that writing the successive positions of a particle on such an 
orbit every time the bars align generates a closed curve, which we call a loop. Loops 
allow us to verify consistency of the potential. As maps of doubly periodic orbits, loops 
can be used to search the phase-space in double bars in order to determine the fraction 
occupied by ordered motions. 

1 Introduction 

Bars within bars appear to be a common phenomenon in galaxies. Recent surveys 
show that up to 30% of early-type barred galaxies contain nested bars The 
relative orientation of the two bars is random, therefore it is likely that the 
bars rotate with different pattern speeds. Inner bars, like large bars, are made 
of relatively old stellar populations, since they remain distinct in near infrared 
0. Galaxies with two independently rotating bars do not conserve the Jacobi 
integral, and it is a complex dynamical task to explain how such systems are 
sustained. To account for their longevity, one has to find sets of particles that 
support the shape of the potential in which they move. Particle motion in a 
potential of double bars belongs to the general problem of motion in a pulsating 
potential [HI 0) of which the restricted elliptical 3-body problem is the best 
known example. Families of closed periodic orbits have been found in this last 
problem, where the test particle moves in the potential of a binary star with 
components on elliptical orbits ,2|. However, such families are parameterized by 
values that also characterize the potential (i.e. ellipticity of the stellar orbit and 
the mass ratio of the stars), and their orbital periods are commensurate with 
the pulsation period of the potential. For a given potential, these families are 
reduced to single orbits separated in phase-space. The solution for double bars 
is formally identical, and there an orbit can close only when the orbital period 
is commensurate with the relative period of the bars. Such orbits are separated 
in phase-space, and therefore families of closed periodic orbits are unlikely to 
provide orbital support for nested bars. Another difficulty in supporting nested 
bars is caused by the piling up of resonances created by each bar, which leads 
to considerable chaotic zones. In order to minimize the number of chaotic zones. 
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resonant coupling between the bars has been proposed pi) , so that the resonance 
generated by one bar overlaps with that caused by the other bar. 

Finding support for nested bars has been hampered by the fact that closed 
periodic orbits are scarce there. However, it is particles, not orbits, which create 
density distributions that support the potential. The concept of closed periodic 
orbit is too limiting in investigation of nested bars, and another description of 
particle motion, which does not have its limitations, is needed. Naturally, in 
systems with two forcing frequencies, double-periodic orbits play a fundamental 
role. Thus in double bars a large fraction of particle trajectories gets trapped 
around a class of double-periodic orbits. Although such orbits do not close in 
any reference frame, they can be conveniently mapped onto the loops which 
are an efficient descriptor of orbital structure in a pulsating potential. The loop 
is a closed curve that is made of particles moving in the potential of a doubly 
barred galaxy, and which pulsates with the relative period of the bars. Orbital 
support for nested bars can be provided by placing particles on the loops. 

Here I give a systematic description of the loop approach, which recovers 
families of stable double-periodic orbits, and which can be applied to any pul- 
sating potential. In §2 I use the epicyclic approximation to introduce the basic 
concepts, and in §3 I outline the general method. 

2 The epicyclic solution for any number of bars 

If a galaxy has a bar that rotates with a constant pattern speed, it is convenient 
to study particle orbits in the reference frame rotating with the bar. If two or 
more bars are present, and each rotates with its own pattern speed, there is 
no reference frame in which the potential remains unchanged. In order to point 
out formal similarities in solutions for one and many bars, I solve the linearized 
equations in the inertial frame. This is equivalent to the solution in any rotating 
frame, and the transformation is particularly simple: in the rotating frame the 
centrifugal and Coriolis terms are equivalent to the Doppler shift of the angular 
velocity. It is convenient to show it in cylindrical coordinates {R^ip,z): if is 
the rotation axis, then the R and (p components of the equation of motion for 
the rotating frame, = — — 2(J7b x r) — x (rJe x r), can be written as 



These equations are identical with the components of the equation of motion in 
the inertial frame. 



Rip + 2R{^ + {2b) 



R-R{p> + QBf 



dR' 
1 

Rd^ 



r = -V^, 



(1) 



where clearly the angular velocity ip in the rotating frame corresponds to ip + {1b 
in the inertial frame. For the rest of this section I assume the inertial frame, in 
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which the equation of motion (^) has the following R and f components in 
cylindrical coordinates 

1 d<P 

^^^'^^-~rW 

The z component in any frame is z = —d<P/dz, but I consider here motions in 
the plane of the disc only, hence I neglect the dependence on z. 

To linearize equations (O and (j^J, one needs expansions of i?, ip and <P to first 
order terms. The epicyclic approximation is valid for particles whose trajectories 
oscillate around circular orbits. For such particles one can write 

R{t) =Ro + Ri{t), (4) 
(p{t) ^ ipQQ + [2ot + ipi{t), (5) 
'P{R,^,t)='Po{R) + ^i{R.V,t), (6) 

where terms with index I are small to the first order, and second- and higher- 
order terms were neglected. The parameter (ySpo allows the particle to start from 
any position angle at time t = 0, so that — ipoo + Q^t. Asymmetry in the 
potential is small and may be time-dependent. The angular velocity f2o on the 
circular orbit of radius i?o relates to the potential through the zeroth order 
of HJ: 12^ = {l/RQ){d<Po/dR) \r„. The zeroth order of © is identically equal to 
zero, and the first order corrections to and Q take respectively forms 

d<pT 

Ri - AA^oRi - 2Rono^i = (7) 

1 d<Pj 

RqlPi + 2QoRi ^^^"9^'^°''^°' 

where A is the Oort constant defined by 4Ai7o = ~ ^q§t \ r,, ■ 

We assume that the bars are point-symmetric with respect to the galaxy 
centre. Thus to first order the departure of the barred potential from axial sym- 
metry can be described by a term cos{2ip). If multiple bars, indexed by i, rotate 
independently as solid bodies with angular velocities f2i, the time-dependent 
first-order correction $j to the potential can be written as 

^iiR,ip,t) = Y,^^iR)cos[2i^ ^ f2,t)], (9) 

i 

where the radial dependence ^i{R) has been separated from the angle depen- 
dence. No phase in the trigonometric functions above means that we define t — Q 
when all the bars are aligned. Derivatives of © enter right-hand sides of {Tj) and 
(|HJ), which after introducing LOi = 2{Q[) — Qi) take the form 

Ri - AAQ^Ri - 2Rof2oipi = - 1] ^ 1% co8{Lo,t + 2</joo), (10) 



2 

Rq 



Roifi + 2QqRi =-5-2^ '^'(^0) sin(^'^ + 2'^5oo). (11) 
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In order to solve the set of equations (|1UI11|I . one can integrate and 
get an expression for Rotpi, which furthermore can be substituted to H10() . This 
substitution eUminates (/?/, and one gets a single second order equation for Rj, 
which can be written schematically as 

Ri + kIRi = ^ cos(w,t + 2(poo) + C^, (12) 

i 

where Ai = -^^^ ~ W|i?o' '^o = 4J7o(i7o - A), and C^/2Qq is the integration 
constant that appears after integrating (|ll(l . This is the equation of a harmonic 
oscillator with multiple forcing terms, whose solution is well known. It can be 
written as 

Ri{t) = Ci cos(Koi + ^) + X! cos(cjit + 2(^oo) + C^/nl. (13) 

i 

The first term of this solution corresponds to a free oscillation at the local 
epicyclic frequency kq, and Ci is unconstrained. The terms under the sum de- 
scribe oscillations resulting from the forcing terms in @, and Mi are functions 
of Ai. Hereafter I focus on solutions without free oscillations, thus I assume that 
Ci = 0. These solutions will lead to closed periodic orbits and to loops. The 
formula for ^piit) can be obtained by substituting (|13|l into the time-integrated 
l(TT|l. As a result, one gets 

2AC 

(^/ = ViV,cos(oj,t + 2v?oo)-^^, (14) 

where again Ni are determined by the coefficients of the equations above. Note 
that to the first order i7o[i?o + C^/n^] = i7o[i?o] — 'i.AC^/ KqRq, thus the inte- 
gration constants entering H13I) and 1141) correspond to a change in the guiding 
radius i?o, and to the appropriate change in the angular velocity Qq. They all 
can be incorporated into i?o, and in effect the unique solutions for i?/ and Lpi 
are 

Ri{t) = ^M, cos(w,;t + 2(^oo), (15) 

i 

ipi{t) = ^ N- cos{uj.it + 2ipoo) + const, (16) 

i 

where free oscillations have been neglected. The integration constant in H16|l is 
an unconstrained parameter of the order of ipi. 

2.1 Closed periodic orbits in a single bar 

In a potential with a single bar there is only one term in the sums (|15ll and 
hereafter indexed with B. Consider the change in values of i?/ and ipj for a 
given particle after half of its period in the frame corotating with the bar. This 
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interval is taken because the bar is bisymmetric, so its forcing is periodic with 
the period tt in angle. After replacing t by i + tt/ {Hq — Qb) one gets 

Ri = Mb cos[ujB{t + -tT^^-Ft) + ^V'oo] 

= Mb cos{LUBt + 2tt + 2</joo). 

Thus the solution for i?/ after time 7r/(f?o — ^?_b) returns its starting value, 
and the same holds true for tpi. After twice that time, i.e. in a full period of this 
particle in the bar frame, the epicycle centre returns to its starting point and the 
orbit closes. Thus (fTK|l and (fT?^ describe closed periodic orbits in the linearized 
problem of a particle motion in a single bar. 



2.2 Loops in double bars 

When two independently rotating bars coexist in a galaxy (hereafter indexed by 
B and S), there is no reference frame in which the potential is constant. Thus 
when a term from one bar in H15(l and (|16(l returns to its starting value, the term 
from the other bar does not (unless the frequencies of the bars are commensu- 
rate). Therefore the particle's trajectory does not close in any reference frame. 
However, consider the change in value of Rj and ipi after time t: /{fis ~ 
which is the relative period of the bars. One gets 

Ri = Mb cos[ws(i + 7^ j^) + 2^oo] + Ms cos[ws(i + 7^ t^) + 2<^oo] 

J/5 — llB S^S — J^S 

= Mb co8{ujBt + ^ + 2</joo) + Ms cos(cJsi + "^t^tT 7^ + 2<^oo) 

Us — Ub J/S — 

= Afs cos(wBt + 27r + 2(pm) + Mscos{ujst + 2(/3oi), 

where (/Jqi = yoo + 7'' f2s-f2B • same result can be obtained for 93/. This means 
that the time transformation t ^ t + 7r/(f?s — f?^) is equivalent to the change 
in the starting position angle of a particle from (/jqo to 9591 ■ Consider motion 
of a set of particles that have the same guiding radius Rq, but start at various 
position angles (^oo- This is a one-parameter set, therefore in the disc plane it is 
represented by a curve, and because of continuity of and (jl6|l this curve is 
closed. After time T:/{f2s — ^b), a particle starting at angle ipoo will take the 
place of the particle which started at ipQi, a particle starting at ipoi will take the 
place of another particle from this curve and so on. The whole curve will regain 
its shape and position every t: /{f2s ^ ^b) time interval, although positions of 
particles on the curve will shift. This curve is the epicyclic approximation to the 
loop: a curve made of particles moving in a given potential, such that the curve 
returns to its original shape and position periodically. In the case of two bars, 
the period is the relative period of the bars, and the loop is made out of particles 
having the same guiding radius Rq. Particles on the loop respond to the forcing 
from the two bars, but they lack any free oscillation. An example of a set of 
loops in a doubly barred galaxy in the epicyclic approximation can be seen in 
Since they occupy a significant part of the disc, one should anticipate large 
zones of ordered motions also in the general, non-linear solution for double bars. 
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3 Full nonlinear solution for loops in nested bars 

Tools and concepts useful in the search for ordered motions in double bars are 
best introduced through the inspection of particle trajectories in such systems. 
For this inspection I chose the potential of Model 1 defined in |H], where the 
small bar is 60% in size of the big bar, and pattern speeds of the bars are not 
commensurate. Consider a particle moving in this potential inside the corotation 
of the small bar. Simple experiments with various initial velocities show that if 
the initial velocity is small enough, the particle usually remains bound. A typical 
trajectory is shown in the left panels of Fig.l - since it depends on the reference 
frame, it is written twice, for reference frame of each bar. Further experimenting 
with initial velocities shows that particle trajectories are often even tidier: they 
look like those in the right panels of Fig.l, as if the trajectories were trapped 
around some regular orbit. 

Fine adjustments of the initial velocity lead to a highly harmonious trajectory 
(Fig. 2), which looks like a loop orbit in a potential of a single bar (see e.g. Fig. 3. 7a 
in PI). This is only a formal similarity, but understanding it will let us find out 
what kind of orbit we see in Fig. 2. The loop orbit in a single bar forms when a 
particle oscillates around a closed periodic orbit. Therefore two frequencies are 
involved: the frequency of the free oscillation, and the forcing frequency of the 
bar. On the other hand, the Fourier transform of the trajectory from Fig. 2 shows 
two sharp peaks at frequencies equal to the forcing frequencies of the two bars 
(Fig. 3). Thus the trajectory from Fig. 2 also has two frequencies: this time these 
are the forcing frequencies from the two bars, while the free oscillation is absent. 
This is how the solution in the linear approximation (§2.2) was constructed. 




Fig. 1. Two example trajectories (one in the two left panels, one in the two right 
ones) of a particle that moves in the potential of two independently rotating bars. The 
particle is followed for 10 relative periods of the bars, and its trajectory is displayed in 
the frame corotating with the big bar (top panels), and the small bar (bottom panels). 
Each bar is outlined in its own reference frame by the dotted line. Large dot marks the 
starting point of the particle. 
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Fig. 2. A doubly periodic orbit in the doubly barred potential, followed for 20 relative 

periods of the bars, and written in the frame corotating with the big bar (left) , and the 
small bar (right). The long axis of each bar is marked by the dashed line. Dots mark 
positions of the particle at every alignment of the bars. 



We conclude that in both the linear (epicyclic approximation) case and the 
general case wc arc dealing with doubly periodic orbits in an oscillating potential 
of a double bar, with frequencies equal to the forcing frequencies of the bars. 
In the epicyclic approximation, these orbits have a nice feature that particles 
following them populate loops: closed curves that return to their original shape 
and position at every alignment of the bars. One may therefore expect that also 
in the general case these particles gather on loops. 

If in the general case particles on doubly periodic orbits form a loop, one can 
construct it by writing positions of a particle on such an orbit every time the 
bars align. These positions are the initial conditions for particles forming the 
loop, because after every alignment, the n*'' particle generated in this way takes 
the position of particle n + 1. The first 20 positions of a particle on a doubly 



freqijericy 

Fig. 3. Fourier transforms of the trajectories from right panels of Fig.l (dotted line) and 
from Fig. 2 (solid line). The peaks in the solid line are related to the forcing frequencies 
of the bars, and the peaks in the dotted line are not. 
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Fig. 4. Evolution of tlie loop from Fig. 2 during one relative period of the bars. The 
bars, outlined with solid lines, rotate counterclockwise. The loop is made out of points 
that represent separate particles on doubly periodic orbits. 

periodic orbit are overplotted in Fig. 2. They indeed seem to be arranged on a 
closed ellipse-like curve; the shape of this curve varies in time, but it returns to 
where it started at every alignment of the bars (Fig. 4). This construction shows 
that in the general case particles on doubly periodic orbits also form loops. Note 
that positions of particles on other orbits, which involve free oscillations, when 
written at every alignment of the bars, densely populate some two-dimensional 
section of the plane, and do not gather on any curve. It is extremely useful for 
the investigation of the orbital structure in double bars that the appearance of 
the loop is frame- independent. Loops provide an efficient way to classify doubly 
periodic orbits, which has been hampered so far by the dependence of the last 
ones on the reference frame. 

It turns out that doubly periodic orbits play crucial role in providing orbital 
support for the pulsating potential of double bars. No closed periodic orbits 
have been proposed as candidates for the backbone of such a potential. If in a 
given potential of two bars there are loops that follow the inner bar, and other 
loops that follow the outer bar, then one may expect that such a potential is 
dynamically possible. An example of such a potential has been constructed in 
|H]. The loop from Fig. 4 does not follow either bar in its motion, and therefore 
it is unlikely that it supports the assumed potential. It can be shown that in 
that potential, there are no loops which could support the two bars. Thus that 
potential is not self-consistent. This example shows how efficient is the loop 
approach in rejecting hypothetical doubly barred systems that have no orbital 
support. 

Doubly periodic orbits in double bars are surrounded by regular orbits in 
the same way as are the closed periodic orbits in a single bar. In both cases, 
the trapped regular orbits oscillate around the parent orbit. The trajectory from 
the right panels of Fig.l is an example of a regular orbit that is trapped around 
the doubly periodic orbit from Fig. 2. How much of the phase space in double 
bars is occupied by orbits trapped around doubly periodic orbits? It can be 
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starting position 



Fig. 5. The width of the ring formed by particles trapped around loops in Model 2 
from [H] as a function of the particle's position along the minor axis of the aligned bars, 
and of its velocity (perpendicular to this axis). Darker color means smaller width. In 
the insert, the same is shown for rings around closed periodic orbits in a single bar 
(same model, but inner bar axisymmetric) . Regions related to the Xi and X2 orbits, 
and to the loops originating from them, are marked. 

examined by launching a particle from e.g. the minor axis of the bar, in the 
direction perpendicular to this axis, when the bars are aligned. If the particle 
is trapped, its positions at every alignment of the two bars lie within a ring 
containing the loop. The width of this ring depends on the particle's position 
along the minor axis, and on its velocity. It is displayed in Fig. 5 for the potential 
of Model 2 defined in |S] . Two stripes of low width appear on the diagram, which 
correspond to the xi and X2 orbits in a single bar (displayed in the insert). 
Thus in double bars there are doubly periodic orbits that correspond to closed 
periodic orbits in single bars. There are possible regions of chaos in double bars 
(white stripes in Fig. 5), but overall loops in double bars and periodic orbits in 
single bars trap similar volumes of phase-space around them. 



4 Conclusions 

In a potential of two independently rotating bars, a large fraction of phase space 
can be occupied by trajectories trapped around parent regular orbits. These 
orbits are doubly periodic, with the two periods corresponding to the forcing 
frequencies of the two bars, but they do not close in any reference frame. Like 
particle trajectories oscillating around closed periodic orbits in a single bar, 
particle trajectories in double bars oscillate around the doubly periodic parent 
orbits. The structure of the parent regular orbits can be mapped using the loop 
approach, which allows us to single out dynamically possible double bars. 
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